Reentrant fractional quantum Hall states in bilayer graphene: Controllable, driven 

phase transitions 
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Here we report from our theoretical studies that in biased bilayer graphene, one can induce phase 
transitions from an incompressible state to a compressible state by tuning the bandgap at a given 
electron density. Likewise, variation of the density with a fixed bandgap results in a transition 
from the FQHE states at lower Landau levels to compressible states at intermediate Landau levels 
and finally to FQHE states at higher Landau levels. This intriguing scenario of tunable phase 
transitions in the fractional quantum Hall states is unique to bilayer graphene and never before 
existed in conventional semiconductor systems. 



The unconventional quantum Hall effect in monolayer 
graphene, whose experimental observation 0] unleashed 
quite unprecedented interest in this system Q, reflects 
the unique behavior of massless Dirac fermions in a mag- 
netic field [J, [|| . In bilayer graphene this effect confirms 
the presence of massive chiral quasiparticles [6] . An im- 
portant characteristic of bilayer graphene is that it is a 
semiconductor with a tunable bandgap between the va- 
lence and conduction bands 7]. This modifies the Lan- 
dau level spectrum and influences the role of long-range 
Coulomb interactions Q. Here we report that the frac- 
tional quantum Hall effect (FQHE), a distinct signature 
of interacting electrons in the system [j| is very sensi- 
tive to the interlayer coupling strength and the bias volt- 
age. We propose that by tuning the bias voltage, one can 
induce phase transitions from an incompressible state to 
a compressible state at a given gate voltage. In the same 
vein, variation of the gate voltage at a fixed bias poten- 
tial results in a transition from the FQHE states in lower 
Landau levels to compressible states in intermediate Lan- 
dau levels and finally to FQHE states in higher Landau 
levels. This interesting scenario of tunable phase tran- 
sitions in the FQH states is unique to bilayer graphene 
and never before existed in conventional semiconductor 
systems. The FQHE in mon olay er graphene was in fact, 
studied theoretically by us [Tl| and subsequent experi- 
ments confirmed the existence of that effect in suspended 
monolayer graphene samples [Hj]. No such studies have 
been reported on bilayer graphene. 

We assume that the bilayer graphene consists of two 
coupled graphene layers with the Bernal stacking ar- 
rangement. In that case, we are mainly concerned with 
the coupling between the atoms of sublattice A of the 
lower layer and atoms of sublattice B' of the upper 
layer. The single-particle levels in bilayer graphene have 
two-fold spin degeneracy and two-fold valley degeneracy, 
which can be lifted in many-particle systems at relatively 
large magnetic fields Considering only one valley 

(say valley K) and one direction of spin, we describe the 
state of the bilayer system in terms of the four-component 



spinor (tpA, ipB, ipB' , *pA') T ■ Here subindices A, B and A', 
B' correspond to lower and upper layers respectively. The 
strength of inter-layer coupling is described in terms of 
the inter-layer hopping integral, t. In a biased bilayer 
graphene the bias potential is introduced as the poten- 
tial difference, At/, between the upper and lower layers. 
The Hamiltonian of the biased bilayer system in a per- 
pendicular magnetic field then takes the following form 
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where ir± = n x ± in. 



= p + eA/c, p is an electron 

two-dimensional momentum, A is the vector potential, 
and v F ~ 10 6 m/s is the fermi velocity. 

In a perpendicular magnetic field the Hamiltonian 
([T]) generates a discrete Landau level energy spectrum. 
The corresponding eigenfunctions can be expressed in 
terms of the conventional nonrelativistic Landau func- 
tions. The electron states at sublattices A and A' are 
expressed in terms of the n-th 'nonrelativistic' Landau 
functions, while the electron states at sublattices B and 
B' are described by the \n — 1| and n + 1 Landau func- 
tions, respectively. Therefore the Landau states in bi- 
layer graphene can be described as a mixture of the n, 
n + 1, and n — 1 nonrelativistic Landau functions be- 
longing to different sublattices Q. This mixture, for a 
given value of n, results in four different Landau levels 
in bilayer graphene. The energies, e, of the four Landau 
levels corresponding to the index n can be found from 
the following equation Q 



(e + 5f - 2(n + 1) [(e - 5f - 2n] = (e 2 - S 2 )t 2 , (2) 



where S = AU/2 and all energies are expressed in units 
of hvp/ls- Here Ib = (h/eB)? is the magnetic length. 

It is convenient to introduce special labeling of the Lan- 
dau levels in bilayer graphene. From Eq. ([2]), we can see 
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that for each value of n, n = 0,1,2,..., there are four 
solutions, i.e., four Landau levels. Then each of these 
Landau levels can be labeled as nu), where i = 1,...,4 
is the number of the Landau level corresponding to the 
solution of Eq. @ for a given value of n. 

For a partially occupied Landau level, the properties 
of the system, e.g., the ground state and excitations, are 
completely determined by the inter-electron interactions, 
which can be expressed by Haldane's pseudopotentials, 
V m , [HI which are the energies of two electrons with rel- 
ative angular momentum m. In a graphene bilayer the 
Haldane pseudopotentials in a Landau level with index 
n and the energy e have the form 



Vt ] =i %V{q)[F n ,M] 2 L m {q 2 )e 



2tt 



(3) 



where L m (x) are the Laguerre polynomials, V(q) — 
2ire 2 1 (kIbQ) is the Coulomb interaction in the momen- 
tum space, k is the dielectric constant, and F n ^ £ (q) are 
the corresponding form factors 
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The form factors of bilayer graphene [Eq. (UJ)] are clearly 
different from the corresponding form factors of a sin- 
gle layer of graphene. For a single layer of graphene, 
the form factor of the n — Landau level is the same 
as that of conventional non-relativistic electrons [ll|> [HI > 
F n —o(q) = Lo (g 2 /2). The form factors of higher Landau 
levels are determined by the mixture of L n and L n+ i 
terms. For bilayer graphene the form factors of the Lan- 
dau level with index n = are the mixtures of the Lq 
and h\ terms and are different from that in the non- 
relativistic case. There is of course, one special Landau 
level in bilayer graphene with index n = 0, whose proper- 
ties are completely identical to that of the non-relativistic 
n = Landau level. Indeed, it is clear from Eq. ([2]) that 
at n = there is a Landau level with energy e = S. This 
energy does not depend on the coupling between the lay- 
ers, i.e., the inter-layer hopping integral, t. The form 
factor of this Landau level is exactly equal to the form 
factor of a non-relativistic system of the n = Landau 
level, F n =o,e=6 = Lq. Therefore, all many-body proper- 
ties of a bilayer system in the n = 0, e = 6 Landau level 
are completely identical to those of a non-relativistic con- 
ventional system in the n = Landau level. 

For Landau levels with higher indices the form fac- 
tor is the mixture of three different functions, L„, L n _i, 
and L n+ \. Therefore, in general, the strength of inter- 
electron interactions in bilayer graphene is strongly mod- 
ified as compared to its value in monolayer graphene. To 



address the effects of these modifications on the proper- 
ties of the many-electron system in a graphene bilayer we 
study below the partially occupied Landau levels with 
fractional filling factor corresponding to the FQHE (9[. 
In what follows, we study the many-electron system at 
various fractional fillin g fa ctors numerically within the 
spherical geometry [Hl,Il3]. The radius of our spehere is 
R = vSIb, where 25* is the number of magnetic fluxes 
through the sphere in units of the flux quanta. The 
single-electron states are characterized by the angular 
momentum, 5, and its z component, S z . For the many- 
electron system the corresponding states are classified 
by the total angular momentum L and its z component, 
while the energy of the state depends only on 
given fractional filling of the Landau level is determined 
by a special relation between the number of electrons N 
and the radius of the sphere R. For example, the 
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FQHE state is realized at S = (|)(iV - 1), while the 



|-FQHE state corresponds to the relation S 
With the Haldane pseudopotentials [Eq. ([2 
mine the interaction Hamiltonian matrix [1 



(l)N-2. 
wc dctcr- 



and then 

calculate a few lowest eigenvalues and eigenvectors of this 
matrix. The FQHE states are obtained when the ground 
state of the system is an incompressible liquid, the energy 
spectrum of which has a finite many-body gap 0, H3] • 

We begin our study with the celebrated i-FQHE [To| 
corresponding to the filling factor v = |. In Fig. [T]we 
show the behavior of the Landau level spectra as a func- 
tion of the bias voltage and for different values of the 
inter-layer hopping integral, t. Only the Landau levels 
with positive energies (corresponding to the conduction 
band) arc shown in the figure. A similar behavior is valid 
for other FQHE filling factors, e.g., for v = ~. Figure Q] 
clearly illustrates that the FQHE can be observed in all 
n = Landau levels with the strongest FQHE being in 
the first and the third n = Landau levels, i.e., O^) and 
0(3). An interesting behavior was also observed for the 
n = 1 Landau levels. For AU w 200 meV the n = 1 Lan- 
dau levels show an anti-crossing, which is accompanied 
by strong changes in the properties of the FQHE. More 
specifically, at a small AU the FQHE can be observed 
only in the lower n = 1 Landau level, while for larger 
AU the FQHE is possible only in higher n = 1 Lan- 
dau level. These have important implications for pos- 
sible experimental observations of this unique behavior 
(see Fig. □»): 

(i) By applying a gate voltage the electron density can 
be tuned so that the first three Landau levels in bilayer 
graphene are completely occupied and the next Landau 
level is partially occupied with the FQHE filling factor, 
for example, v = According to the notations of Fig. [TJ 
this means that the 0(i), 0(2), 0(3) Landau levels are fully 
occupied, while the In) Landau level has a filling factor 
■|. Then, by varying AU from a small value, e.g., 100 
meV, to a larger value, e.g., 400 meV, one can observe 
the disappearance of the FQHE (see line (i) in Fig. [Jib)). 

(ii) The bias voltage is kept fixed at a large value, e.g., 
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FIG. 1: A few lowest Landau levels of the conduction band 
are shown as a function of the bias potential, AU, for dif- 
ferent values of inter-layer coupling: (a) t = 100 meV (b) 
t — 200 meV and (c) t = 400 meV and the magnetic field is 
15 Tesla. The numbers next to the curves denote the corre- 
sponding Landau levels. The Landau levels where the FQHE 
can be observed are drawn as blue and green dots. The green 
dots correspond to the Landau levels where the FQHE states 
are identical to that of a monolayer of graphene or a non- 
relativistic conventional system. The red dots represent Lan- 
dau levels with weak FQHE. The open dots correspond to 
the Landau levels where the FQHE cannot be observed. In 
(b), the dashed lines labeled by (i) and (ii) illustrate two sit- 
uations: (i) under a constant gate voltage and variable bias 
potential; (ii) under a constant bias potential and variable 
gate voltage. 



AU = 400 meV. Then by varying the gate voltage and 
thus increasing the electron density, one can observe the 
disappearance and reappearance of the FQHE at large 
Landau levels (see line (ii) in Fig. Q3b)). 

The appearance or suppression of the FQHE in differ- 
ent Landau levels is correlated with the unique behavior 
of Haldane pseudopotentials in bilayer graphene. The 
FQHE is observed in Landau levels with a rapid decrease 
of the corresponding pseudopotentials with increasing an- 
gular momentum, m. For the first and second n = 
Landau levels (0(i) and 0( 2 ) levels), between m = 1 and 




Angular momentum (L) Angular momentum (L) 



FIG. 2: Low energy excitation spectra of the FQHE states 
shown for different values of the bias potential. The numbers 
next to the lines are the values of the bias potential in meV. 
(a) v = i-FQHE (eight electrons) system in the Landau level 
1(2), (b) v = |-FQHE (eight electrons) system in the Lan- 
dau level 1(!), (c) v — |-FQHE (ten electrons) system in the 
Landau level 1(2), (d) v = |-FQHE (ten electrons) system in 
the Landau level lm. The systems are fully spin-polarized. 
The flux quanta is 25 = 21. The magnetic field is 15 T. The 
ground state is shown as a solid dot at L = 0. 

to = 3, the pseudopotential has the fastest decay for the 
first level, which results in a higher FQHE gap in this 
Landau level. The pseudopotentials in the 0(2) Landau 
level exactly coincide with the pseudopotentials of the 
?7 = Landau level in a single layer of graphene. For 
the 77 = 1 Landau levels, the pseudopotentials in one of 
the 77 = 1 levels are close to those in the n = 1 Landau 
level of the monolayer graphene, resulting in a well pro- 
nounced FQHE. While in the other n = 1 Landau level 
of the bilayer graphene the pseudopotentials display a 
slow decrease with m and no FQHE is expected in that 
Landau level. 

The collapse of the FQHE gap, corresponding to the 
appearence of anticrossing of the n — 1 Landau levels, 
is illustrated in Fig. [3] The FQHE gap has a monotonic 
dependence on the bias voltage. In the anticrossing re- 
gion the gap disappears for the lower n = 1 Landau level 
(see Fig.[2k,c) and reappears for the higher 77 = 1 Landau 
level (see Fig. [5b, d). The evolutions of the energy spectra 
of the incompressible liquid with different filling factors 
are similar, which is illustrated in Fig. &p,,b (for v = |) 
and Fig.f2b,d (for v = |). This behavior was never before 
observed in the FQHE of conventional two-dimensional 
electron systems. 

The strength of the FQHE effect, i.e., the magnitude 
of the excitation gap, depends on the parameters of the 
bilayer graphene, i.e., on the bias voltage and the inter- 
layer hopping integral, t. In Fig. [3] such dependence is 
shown for |- and |-FQHE in different Landau levels as 
a function of the hopping integral, t. In accordance with 
the properties of Haldane pseudopotentials, the excita- 
tion gap of the 0(2) Landau levels does not depend on 
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FIG. 3: The FQHE gaps are shown for different Landau lev- 
els. The labels next to the lines correspond to the labeling 
of Landau levels shown in Fig. [1] (a) v = i-FQHE (eight 
electron) system at AU = 10 meV, (b) v = i-FQHE (eight 
electron) system at AU = 400 meV, (c) v = |-FQHE (ten 
electron) system at AU = 10 meV, (d) v = |-FQHE (ten 
electron) system at AU = 400 meV. All systems are fully 
spin polarized. The flux quanta is 2S = 21. The magnetic 
field is 15 Tesla. 



the bias voltage and on the inter-layer hopping integral. 
The corresponding gap remains constant and is equal to 
gap of the FQHE in a single layer of graphene in the n = 
Landau level. This gap has the smallest value compared 
to that of the other Landau levels. For zero hopping inte- 
gral the two layers of graphene become decoupled and the 
bilayer system becomes identical to a single layer with ad- 
ditional double degeneracy. This property is clearly seen 
in Fig. [31 where at t = there are only two doubly de- 



generate FQHE gaps, corresponding to n — and n = 1 
single layer Landau levels. 

At a small bias voltage (AU = 10 meV, see Fig. [3^,c), 
the FQHE gaps have a monotonic dependence on the 
hopping integral. The values of the gaps increases for 
the n = Landau levels, while for the n = 1 Landau 
level the gap decreases. At large values of the bias volt- 
age, AU = 400 meV (see Fig. [3b,d), the FQHE gap in 
the O3 Landau level has a nonmonotonic dependence on 
the hopping integral with a maximum around i = 100 
meV. These results also illustrate that the FQHE in bi- 
layer graphene can be more stable, i.e., the corresponding 
FQHE gap is much larger than in the case of FQHE in a 
single graphene layer. As an illustration, the FQHE gaps 
in the O^j and 0( 3 ) Landau levels at finite values of t are 
larger than the FQHE gaps in a single layer of graphene, 
i.e., at t = 0. 

Our present work suggests an unique opportunity to 
tune the reentrant FQHE states. The transition from 
the FQHE state to a gapless compressible state is con- 
tinuous. The gap decreases monotonically with the bias 
potential and finally collapses. The collapse of the FQHE 
gap occurs at a non-zero angular momentum, L rs 2, re- 
sulting in the formation of the ground state with a finite 
momentum. Experimental observation of this transition 
will provide unique opportunities, for the first time, to 
study the phase transition from an incompressible liq- 
uid to a possible charge density wave state and then to 
another incompressible state. 
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